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Abstract 

We show that gravitational radiation is produced quite efficiently in interac¬ 
tions of classical waves created by resonant decay of a coherently oscillating 
field. For simple models of chaotic inflation in which the inflaton interacts 
with another scalar field, we find that today’s ratio of energy density in grav¬ 
itational waves per octave to the critical density of the universe can be as 
large as at the maximal wavelength of order 10^ cm. In the pure Xif'^ 

model, the maximal today’s wavelength of gravitational waves produced by 
this mechanism is of order 10® cm, close to the upper bound of operational 
LIGO and TIGA frequencies. The energy density of waves in this model, 
though, is likely to be well below the sensitivity of LIGO or TIGA at such 
frequencies. We discuss possibility that in other inflationary models interac¬ 
tion of classical waves can lead to an even stronger gravitational radiation 
background. 
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I. INTRODUCTION 


Recent research in inflationary cosmology has attracted attention to highly non- 
eqnilibrinm states created in a decay of a coherently oscillating held after the end of inflation. 
These states could support a number of non-equilibrium phenomena, such as non-thermal 
symmetry restoration and baryogenesis |^,0] shortly after or during the decay of the 
oscillating held. 

In this paper we want to show that non-equilibrium states produced by the decay of 
coherent oscillations of a held are a quite efficient source of a stochastic background of 
gravitational waves. 

There are several possible processes in the early universe capable to produce a stochastic 
background of relic gravitational waves. One is the parametric amplification of vacuum 
graviton fluctuations during inflation |Q. This process is efficient on all frequency scales. 
Waves with lowest frequencies cause inhomogeneities of cosmic microwave background [Q. 
In conventional scenarios, this restricts the amplitude of high-frequency gravitational waves 
to be far below [Q the experimental limits accessible for direct detection experiments in near 
future ; this conclusion changes in superstring motivated cosmologies 0. Another source of 
gravitational radiation is classical emission that accompanies collisions of massive bodies. A 
natural source in this class in the early universe is a strongly first-order phase transition when 
gravitational waves are produced in collisions of bubbles of a new phase M, in particular 
the phase transition that terminates first-order inflation |^. Gravitational radiation is also 
emitted during the decay of a cosmic string network 0- 

In this paper we discuss a new source of relic gravitational waves. Decay of coherent os¬ 
cillations of a scalar field can produce large, essentially classical fluctuations via parametric 
amplification (resonance) H- These classical fluctuations, which can be viewed either as 
classical waves traveling through the universe or, at least qualitatively, as quantum “parti¬ 
cles” in states with large occupation numbers, interact with the oscillating background and 


each other. This interaction, which we call rescattering , is accompanied by gravitational 
radiation. That is the effect we want to estimate. 

Favorable conditions for an effective parametric resonance in cosmology naturally appear 
in inflationary models |1I^,pA|. We consider two types of simple inflationary models here. 


One type of models has two scalar fields with an interaction potential of the form g'^cfP'X"^/2, 
and the resonance produces mostly fluctuations of a scalar field X other than the field 0 that 
oscillates (although subsequent rescattering processes produce large fluctuations of the field 
0 as well). We consider a range of moderate values of the coupling (see below); in this 
case fluctuations are not suppressed too strongly by non-linear effects (cf. Refs. jlb HI^] )- In 
simplest models of this type, 0 is the inflaton itself. For these models, we find that typically 


2 













~ 10“® of the total energy of the universe go into gravitational waves, at the time of their 
production. The minimal today’s frequency /min of these waves is typically of order 10^ Hz, 
and today’s spectral density at this frequency can be as large as 10“^^ of the critical density. 

Another type of model we considered was the pure model of chaotic inflation. In 
this model, the minimal today’s frequency /min is of order 10^ Hz, close to the upper bound 
of the operational LIGO and TIGA frequencies |^: 10 Hz < /ligo ~ 10^ Hz. We do not 
yet have efficient means of extrapolating our numerical results for today’s spectral intensity 
to the minimal frequency of this model, but we do not expect it to be above 10“^^ of the 
critical density. That would be well below the sensitivity of LIGO or TIGA at frequencies 
of order 10^ Hz. 

We believe, however, that in the absence of a commonly accepted specihc inflationary 
model, it is premature to rule out possibility of experimental detection, even already by 
LIGO or TIGA, of gravity waves produced by the mechanism we consider here. At the 
end of this paper, we discuss possibility of a stronger background of gravitational waves in 
models with more helds or more complicated potentials. 


II. GENERAL FEATURES OF POST-INFLATIONARY DYNAMICS 


In scenarios where inflation ends by resonant decay of coherent held oscillations, post¬ 
inhationary dynamics has two rapid stages. At the hrst of these, called preheating 


huctuations of Bose helds interacting with the oscillating held grow exponentially fast, as a 
result of parametric resonance, and achieve large occupation numbers. At the second stage, 
called semiclassical thermalization |^, rescattering of produced huctuations smears out the 
resonance peaks in power spectra and leads to a slowly evolving state, in which the power 
spectra are smooth [^,^,^]. The system begins to exhibit chaotic behavior characteristic 
of a classical non-linear system with many degrees of freedom. In the course of subsequent 
slow evolution, the power spectra propagate to larger momenta; we expect that eventually 
this process will lead to a fully thermalized state (which does not admit a semiclassical 
description). 

In many models, we hnd that during the stage of semiclassical thermalization, or chao- 
tization, huctuations grow somewhat beyond their values at the end of the resonance stage. 
In such cases, the most ehective graviton production takes place at the end of and shortly 
after the chaotization stage. 

In order not to conhne ourselves to any particular type of inhationary scenario, we will 
not assume that the oscillating held / is the inhaton itself. (It is distinct from the inhaton, for 
example, in hybrid inhationary models |]^.) Let us denote the amplitude of the oscillating 
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field 0 at the end of the chaotization stage as 0c/i, the frequency of its oscillations as m, and 
the Hubble parameter at that time as Hch- These same parameters at the end of inflation, 
when the oscillations start, will be denoted as 0(0), m(0), and -H(O). In simplest models 
of chaotic inflation, in which 0 is the inflaton, 0(0) ~ Mpi, and -ff(O) ~ m(0). Although 
we will use such models for illustrative purposes, our general formulas do not assume these 
conditions. Oscillations of 0 cannot start unless m(0) ^ hf(0), but on the other hand, 
they can start at m(0) ^ hf(0), if they have to be triggered by some other held. Even if 
hf(0) ~ m(0), as in simple models of chaotic inhation, we still have 

Hch < Hr <^m ( 1 ) 


where Hr is the Hubble parameter at the end of the resonance stage. This is because the 
frequency of oscillations redshifts slower (if at all) than the Hubble parameter. We will use 
the condition (|^ in what follows. 

We will consider models, in which oscillating held 0 interacts with a massless scalar held 
X. The Lagrangians for these models are of the form 

L = ^^'^"0^00.0 + ^g^’^d^Xd.X - H(0, X) (2) 

where 


1/(0,X) = 14(0) + ^/02x2 (3) 

and 10, is a potential for the held 0. We will consider two types of 10: 10 = (massive 

0) and 10 = |A0^ (massless 0). The ehect we discuss is not limited to these particular 
models and exists in a wide variety of inhationary models that have a preheating stage. 
Note that in the model with massless 0, the frequency m of oscillations of 0o at the end of 
chaotization is m a/A0 

ch' 

The oscillating 0 amplihes huctuations of X via parametric resonance. An important 
parameter in the problem is the resonance parameter q. Depending on the choice of the 
potential for 0 (see above), q = 5f^0^(O)/4m^ (massive 0) or g = jHK (massless 0). 

In the case of massive 0, it is useful to introduce, in addition to the resonance parameter 
g, the redshifted resonance parameter at the end of the resonance stage 


qr = q 


0^(A) 

02 ( 0 ) 


(4) 


where tr is the time corresponding to the end of the resonance stage, and (j){tr) is the 
amplitude of oscillations of 0o at that time. In the model with massive 0 and massless X, 


parametric resonance can fully develop in an expanding universe if g^ ^ 1 IJl^. Similarly, 
we can introduce 
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(5) 


qch = 


q^lh 

02(0) 


For uniformity of notation, we will sometimes use or Qch instead of q in the case of massless 
0; there is no difference between the three of these in that case. 

Resonant production is most effective for fluctuations of helds that couple to 0 not too 
weakly but also not too strongly, those with ~ 1- For g^ 3> 1, the maximal size of X 
fluctuations is signihcantly suppressed by non-linear effects |T6|-[l8|. Because, for instance. 


superstring models predict a plethora of scalar helds, we expect that some of those will 
have couplings in the optimal range. So, in what follows we consider moderate values of g,., 
l^qr< 100 . 

In the model with massive 0, the oscillating zero-momentum mode 0o drops rapidly at 
the end of the chaotization stage, and all its energy at that time is transferred to huctuations 


TTI . The variance of 0, ((50)^) at the end of chaotization is thus much larger than 02^. This 


does not happen in the model with massless 0. So, in what follows we consider two cases: 

((50)^) ^ (^Ih and {{6(j)f) > 


For g > 1 in the model where a massless 0 interacts with X, parametric resonance for 
0 itself is insignihcant. To study a case different in that respect, we consider also the pure 
A0^/4 model, in which 0 decays solely due to self-coupling (g2 = 0 in Eq. (H)). 


III. CALCULATIONAL PROCEDURES 


Because the time scale of processes that give rise to gravitational radiation is much 
smaller than the time scale of the expansion of the universe, the energy of gravitational 
waves can be approximately computed starting from the well-known formulas for flat space- 
time. Total energy of gravitational waves radiated in direction n = k/o; in flat space-time 
is [go 

dE 


dfl 


= 2GA,,,im{n) / a;"r^*(k,a;)T'™(k,o^)da; 


( 6 ) 


where T^-^^k, a;) are Fourier components of the stress tensor, and is a projection tensor 
made of the components of n and Kronecker’s deltas. 

For models with two helds (0 and X) with potentials g]) and moderate qr, we will 
present both analytical estimates and numerical calculations based on Eq. (|^). For analyti¬ 
cal estimates, we choose a specihc process—gravitational bremsstrahlung that accompanies 
creation and annihilation of huctuations of the held 0. This is not the only process that can 
produce a signihcant amount of gravitational radiation. For example, because the mass of 
X huctuations oscillates due to interaction with 0o, the collection of X huctuations works 
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(until 00 decays completely) as a gravitational antenna. The reason why we concentrate on 
bremsstrahlung from 0 is that it is a signihcant source of gravity waves with small frequen¬ 
cies, which have the best chance to be detected today. Indeed, we will see that both the 
frequency dependence and the overall magnitude of the effect are reasonably close to those 
of the full intensity for small frequencies that we obtain numerically. So, bremsstrahlung 
from 0 appears to be at least one of the main sources of gravitational radiation with small 
frequencies in the states we consider here. 

Numerical calculations were done for simple models of chaotic inflation, in which 0 is 
the inflaton itself. The calculations were done in conformal time, in which the evolution 
of the system is Hamiltonian (for more detail, see Ref. JT^), and then rescaled back to the 
physical time. We computed directly Fourier transforms cn) over successive intervals 

of conformal time, the duration of which was taken small compared to the time scale of 
the expansion yet large enough to accommodate the minimal frequency of gravity waves 
we had in these calculations. Intensities from different intervals were summed up with the 
weight that takes into account expansion of the universe (see more on that below). Being 
a finite size effect, the minimal frequency of our numerical calculations was much larger 
than the actual minimal frequency of gravity waves, cUmin ~ Hch- hi the models with two 
helds (0 and X), the spectrum obtained numerically is approximated reasonably well by the 
bremsstrahlung spectrum at small frequencies. This allows ns to extrapolate the numerical 
results to a; ~ cUniin using the frequency dependence of bremsstrahlung. 

For the pure A0'^/4 model, the analytical method that we use to estimate the intensity 
of radiation does not apply at frequencies for which we have numerical data. So, we do 
not have efficient means of extrapolating our numerical results to cUmin. For this model, we 
contented ourselves with numerical simulations. 

IV. ANALYTICAL ESTIMATES 

In general, X and 0 fluctuations produced by parametric resonance scatter off the homo¬ 
geneous oscillating background (condensate) of 0, knocking 0 out of the condensate and into 
modes with non-zero momenta. In cases when resonance amplihes mostly fluctuations of X, 
the scattering process can be viewed as a two-body decay of X, X]^ X^' + 0p, in the 

time-dependent background held of the condensate. It can also be thought of as evaporation 
of the inhaton condensate. There is also the inverse process (condensation). 

As we noted in the introduction, one can describe the states we are considering either 
as collections of interacting classical waves, or as collections of “particles” in modes with 
large occupation numbers. For estimating the intensity of bremsstrahlung that accompanies 
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creation and annihilation of fluctuations of the held 0, the description of these huctuations 
in terms of particles is more convenient. For the notion of particle, i.e an entity moving 
freely between collisions, to have anything but purely qualitative meaning, the energy of a 
freely moving “particle” should not be modulated too strongly by its interactions with the 
background. 

In the models (H) with ~ 1, huctuations of 0 are reasonably well described as particles, 
even in the case of massless 0 where they are coupled to the oscillating 0o. Indeed, typical 
momenta p of 0 huctuations at the end of the chaotization stage are of order m, the frequency 
of oscillations of 0o at that time. The frequency of oscillations of a mode with momentum 
p ~ m (the would-be energy of a particle) is only moderately modulated by its coupling to 
the oscillating 0o. Note also that for ~ 1, fhe Hartree correction to the frequency squared 
of 0, does not exceed itself [O, so interaction with X also does not modulate 


frequencies of huctuations of 0 too much. 

So, let us consider huctuations of 0 as particles, neglecting modulations of their energy, 
and estimate the intensity of gravitational bremsstrahlung emitted by these particles in the 
scattering processes described above. We will assume that bremsstrahlung from 0 is at least 
one of the main sources of gravity waves with small frequencies (we will specify the notion 
of “small” frequency below). For estimation purposes, we will neglect its interference with 
possible small-frequency radiation from other sources. 

To describe radiation with small enough frequencies, we can regard the particles—quanta 
of 0—as semiclassical wave packets. This can be seen by considering the Fourier transform 
of the stress tensor of a single classical particle. Using identity transformations, we write 
that stress tensor as 


a;T(-^(k, a;) = —: 


pipj Q 


Aiot—i\<iY 


2m 7-00 p°(l — nv) dt 


dt 


(7) 


where p^, r, and v are the particle’s time-dependent four-momentum, position, and velocity. 
Uncertainty that we can tolerate in the position of a particle, to be able to use this formula, is 
of order 1/k. Uncertainty we can tolerate in momentum is of order of the momentum spread 
among the particles, Ap. Thus, the semiclassical approximation applies when uj = k Ap. 
In the models with two helds, at the end of chaotization, Ap ~ m, so Eq. (^ applies for 
a; -C m. 

On the other hand, in the pure A0^/4 model, fluctuations do not grow during chao¬ 
tization, so for the most efficient graviton production we should consider the end of the 
resonance stage. The power spectrum of fluctuations at the end of resonance, and long into 
chaotization, is concentrated in rather narrow peaks, so that Ap -C m. Then, Eq. ap¬ 
plies only at quite low frequencies, much lower than those we had in numerical simulations. 
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As a result, for this model, we could not check our analytical estimates against numerical 
simulations. 

In the rest of this section we consider the models with two helds. For a; -C m, not only 
we can use the semiclassical expression for the stress tensor of 0, but in addition each act of 
evaporation (or condensation) of a 0 particle can be regarded as instantaneous. So, we can 
use the small frequency approximation, familiar from similar problems in electrodynamics 
It amounts to integrating by parts in (|^) and then replacing the exponential with unity. 
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After that, the integral is trivially taken and depends only on the particle’s hnal (or initial) 
momentum. 

To obtain the power V radiated by a unit volume, we substitute the small-frequency 
limit of Eq. (|^) into Eq. (^), multiply the result by the rate at which 0 fluctuations are 
created (or destroyed), and integrate with the occupation numbers n^{p) of 0 fluctuations. 
The rates of the creation and annihilation processes are almost equal (see below). With 
both processes taken into account, we obtain 


dV 

dui 


4G / n^{p)Tl{p)F{v 


p'^dp 

(^ 


( 8 ) 


where /^(p) is the power spectrum of 0, 7l{p) is the evaporation (condensation) rate, v = 
p/p°, and the function 


F{v) 


\ 3 V 1 — V J 


(9) 


arises after integration over direction of particle’s momentum; we have assumed that both 
the power spectrum and the rate depend only on the absolute value of momentum. Notice 
that the spectrum (Q) is a; independent, as characteristic of small-frequency bremsstrahlung. 

To estimate the rate TZ, we will consider, qualitatively, X fluctuations also as “particles” 
characterized by their own occupation numbers nx{k). This is not a strictly defensible view, 
as X fluctuations are strongly coupled to the oscillating background, but it should do for 
estimation purposes. 

Consider Erst the case ((50)^) ^ The rate Ti{p) can then be estimated as 


p/ N 9^4'lh f d^k nx(k)nx{k-p) 

^ m J (27r)3 Scukcuk.ppo 


( 10 ) 


where a factor of 1/m appears instead of the usual energy delta function because energy of 
the “particles” participating in the process is not conserved, due to time-dependence of 0o. 
This factor estimates the time scale at which energy non-conservation sets in; in our case, 
it is the period of the oscillations of 0o. 










Notice that the rate (p!0[) is larger than the net kinetic rate, which would enter a kinetic 
equation for cj), by a factor of order of a typical occupation number of X. The net kinetic rate 
is the difference between two (large) numbers—the rate at which collisions supply particles 
to a given mode and the rate at which they remove them. But each collision is accompanied 
by bremsstrahlung, so TZ is not the net kinetic rate but a rate of the “in” and “out” processes 
separately. 

For moderate Qr, we can use the following estimates (cf. Refs. |1I^ , |l6 i pT|) : p k ^ ujk ^ 
m, ~ nx ~ 1/5'^- Using these estimates, we obtain lZ{j)) ~ for p m, and 


du 


rr^ 

9W/1, 


Mpi 


( 11 ) 


where in the last relation we have used ~ /(plh dropped a factor of 1/qch- 

An estimate for the total energy density of gravitational waves is obtained by multiplying 
the power V by the time At during which the radiation was substantial. For the ratio of 
Pgw per octave to the total energy density of the universe after chaotization, we obtain 






^Ptot din a;/ 


a; At 


ch 


Mpi 


ujAt 


( 12 ) 


lUpjPtot 

' m‘^(plh ■ The time during which the radiation was 




where in the last relation we used ptot 
substantial is determined by the redshift. At 

For the opposite case {(d0)^) S> we need to replace 0^;^ in the estimate ([l0|) by 
((50)^). Now using ((50)^) ~ and ptot ~ ’^^((d0)^), we see that we need to make the 

same replacement in Eq.(|^. 

These estimates apply at the end of the chaotization stage and extend into the infrared 
up to frequencies of the order of the horizon scale at that time, uj ~ Hch- (We will present 
estimates for today’s pcw below.) 

As an example, consider the model with massive 0 in which 0 is the inflaton itself, and 
m = 10“®Mpi. Take ~ 1, which in this case corresponds to g ~ 10^ |M. The in¬ 


flaton completely decays into fluctuations during the chaotization stage, and at the end 
of it ((50)^) ~ 10“®Mpi 1^. Using ((50)^) in place of in Eq. (p^), we obtain 
PtotdpGw/dlno; ~ 10“® at a; ~ Hch- 

Our analytical estimates are rather crude, as we were somewhat cavalier with numerical 
factors. They do bring out, however, not only the frequency dependence of the effect but to 
some extent also its overall magnitude. Notice, for example, that the last estimate in Eq. 
(0) does not contain powers of or (these two can be expressed through each other 
and cfyl, using the condition qr ~ 1). These quantities, if present, would change the estimate 
by many orders of magnitude, in drastic disagreement with our numerical results. 
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V. BACKGROUND GRAVITATIONAL RADIATION TODAY 


Let us now translate the above estimates into estimates for gravitational background 
radiation today. It is important that gravitational radiation produced after preheating has 
not interacted with matter since then ||22|| . 

First, let us estimate the physical wavelength today, I = 27r/fco, that corresponds 
to a wave vector kch at the end of the chaotization stage. We have fco = Khdch/cbo, 
where Oq is the scale factor today. Thermal equilibrium was established at some tem¬ 
perature T*, and the universe was radiation dominated at that time. This happened 
when relevant reaction rates became equal to the expansion rate. We can find the Hub¬ 
ble parameter at that time as = SnGp^/S, where p* = g^,7i‘^T^/30, and p* = g{T^)] 
g{T) stands for the effective number of ultrarelativistic degrees of freedom at temper¬ 
ature T. The expansion factor from T = T* down to T = Tq ~ 3K is given by 
a*/ao = {go/g^f/^{To/T^) = [g^^/gl^^‘^){3'K^ where go = g{To). De¬ 
pending upon model parameters, the universe could expand from T = T^h to T = T* 
as matter or radiation dominated. We can write = Hchicich/o.*)^^, where a = 2 for 
the radiation dominated case and a = 3/2 for a matter dominated universe. We obtain 
ko = kchach/ao = Khiach/a*){a*/ao) ~ The model de¬ 

pendent factor (ach/a*) does not enter this relation in the radiation dominated case and en¬ 
ters in power 1/4 for a matter dominated universe, i.e. this factor gives a not very important 
correction (at most an order of magnitude). For today’s wavelength, corresponding to wave 
vector kch at the end of chaotization, we hnd I = 27r//co ~ 0.5 {MpiHchy^‘^k~i/{a^:/achY~°‘^‘^ 
cm. 


The smallest wave vector of radiation that could be produced at the end of chaotization 
is of order Hch- The corresponding maximal today’s wavelength /max will fall into the range 
of the LIGO detector, i.e. /max > 3 x 10® cm, when Hch < 10®-10® GeV. 

For illustration, let us consider simple models of chaotic inflation with potentials (^, 
where the oscillating field 0 is the inflaton itself. The Hubble parameter at the end of the 
chaotization stage in these models can be extracted from numerical integrations of Refs. 
r3| , pR|jr7| . For example, for massless inflaton we get Hch/Mpi = (3A/27r)^/^r//^, where Tch 
is conformal time at the end of the chaotization stage. For A = 10“^® and, say, g = 30 we 
get (see Fig. 1) Tch ~ 125. This gives Hch/Mpi ^ 1.4 x 10“^^ and /max ~ 1-3 x 10® cm. 
Another example is the case of massive inflaton with m = 10“®Mpi and q = 10^. Here we 
find Hch/Mpi Ri 2 x 10“®. This gives /max in the range 10^-10® cm. 

An interesting case is that of the pure X(j)‘^/A model (g^ = 0 in (|^)). In this case, 
fluctuations do not grow during chaotization (they actually decrease due to redshift), so 
the most efficient production of gravity waves takes place after the end of the resonance 


10 









(preheating) stage. So, instead of Hch in the above formulas we use Hr, the Hubble parameter 
at the end of resonance. We have Hr/Mp\ = (3A/27r)^/^r“^, where is conformal time at 
the end of resonance. In this model, resonance develops slower that in models with ~ 1- 
As a result. Hr is smaller, and today’s maximal wavelength is larger. For A = 10“^^, we get 
Tr ~ 500 [0, and Hr/Mpi Ri 9 x 10“^^, which gives, for radiation produced at that time, 
^max ~ 5 X 10^ cm. Note, that by the time r 1000 the inflaton still has not decayed in 
this model, fluctuations are still highly non-equilibrium and gravity waves generation will 
continue even at later time. This will increase / ma x by at least another factor of two bringing 
it close to the upper boundary of operational LIGO and TIGA frequencies. 

Now, let us estimate today’s intensity of radiation. The today’s ratio of energy in gravity 
waves to that in radiation is related to (pow/Ptot)c/i via 


Pgw \ _ / Pgw \ / 

Prad / Q \ Ptot \ 9 * / 


(13) 


For the often used parameter ^lg{uj) = {p^rlt^Pcw/d\yiuj)o, where pcrit is the critical energy 
density, we obtain, using ( 0 . 




rsj 


Grad h 


(Plh ^ 



(14) 


where Grad is today’s value of the ratio Prad/Pcrit: Gradh^ = 4.31 x 10 ® When the 
oscillating zero-momentum mode decays completely during the stage of chaotization, in 
(H) is replaced by ((50)^)ch- 

For example, in the case of massless inflaton with A = 10“^^ and q = 30, we have 
10“^Mpj, and with p*/po ~ 100 the estimate (|^) gives G^/i^ ~ 10“^° at / max ~ 10^ 


'Teh 

cm. 


VI. NUMERICAL RESULTS 


We have studied numerically gravitational radiation in the model where p is the massless 
inflaton with = A(/)‘^/4, A = 10“^^ , interacting with a massless held X according to Eq. 
(|). We have used the fully non-linear method of Ref. We have solved equations of 
motion on a lattice 128^ in a box with periodic boundary conditions. 

The model is classically conformally invariant, so by going to conformal variables, the 
equations of motion can be reduced to those in hat space-time. 

Time evolution of the variances (A^) and ((50)^) (in rescaled untis) for resonance pa¬ 
rameter g = 30 is shown in Fig. 1. The rescaled conformal time r is related to time t by 
'/X(j){0)dt = a(r)c/r/a(0). The rescaled conformal helds y and <p are related to the original 
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fields by X = x0(O)a(O)/a(r) and = 00(O)a(O)/a(r). In this model, 0(0) ~ O.SSMpi and 
a(r)/a(0) ~ O.Slr + 1 ||I6| . 

We see that for the above values of parameters, parametric resonance ends at r 73. The 
resonance stage is followed by a plateau (cf. Refs. At the plateau, the variances of 

fluctuations do not grow, but an important restructuring of the power spectrum of X takes 
place. The power spectrum of X changes from being dominated by a resonance peak at 
some non-zero momentum to being dominated by a peak near zero. (For some q, though, 
the strongest peak is close to zero already during the resonance.) When the peak near zero 
becomes strong enough, the growth of variances resumes (in Fig. 1, that happens at r ~ 84), 
and the system enters the chaotization stage. 

The power spectrum of 0 during the time interval of Fig. 1 is shown in Fig. 2. 
The rescaled comoving momentum k is related to physical momentum /Cphys as fcphys = 
A/A0(O)a(O)A;/a(r). Note that in the range 1 ^ /c ^ 10, the power spectrum is approxi¬ 
mately a power law, as characteristic of Kolmogorov spectra p4 . 


Interaction of the helds with gravity is not conformal, and the flat-space-time formula 
(^) for the energy of gravity waves can be used only after we make an approximation. The 
approximation replaces the actual expanding universe with a sequence of static universes. 
Specihcally, conformal time was divided into steps of At = 27tL, where L is the size of the 
integration box, and at each step the physical variables (helds, frequency, and momenta) 
were obtained from the conformal ones, using for a(r) the actual scale factor taken at the 
middle point of the step. The energy of gravitational waves was computed at each step 
using (^, with the corresponding physical variables. Then, energies from all the steps were 
summed up. 

Today’s fig, in physical units, that had been accumulated by conformal time r = 200 
is shown in Fig. 3 for two values of g: q = 30 (dashed line) and q = 105 (dotted line). 
The solid line corresponds to the pure A0^ model (no interaction with X held) and includes 
contributions from times up to r = 1600. In the latter case, the peaks in fig at frequencies 
/ > 4 X 10^ Hz, seen in the hgure, are correlated with peaks in the power spectrum of 0 at 
early stages of chaotization, cf. Ref. [|^|. Since we included contributions to fig from times 
long after the time when the peaks in the power spectrum disappeared, and yet the peaks 
in fig had not been washed out, we suggest that in this model the peaks in fig are a feature 
potentially observable at present. Observation of such peaks could select a particular model 
of inhation. 

We see also that in the model with two interacting helds (0 and X), the linear ut depen¬ 
dence at small uj, characteristic of bremsstrahlung, is overall well born out, both for g = 30 
and g = 105. The minimal today’s frequency in this model is /mm ~ 10® Hz. Extrapolating 
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the results of Fig. 3 to that minimal frequency using the linear law, we obtain the magnitude 
of of order 10“^^ for g = 30 and of order 10“^^ for q = 105. Recall that the analytical 
estimate ([I^) (using ~ 10“^Mpj) gave VLgh? of order 10“^° at the minimal frequency. 

In the pure model, we could not discern any simple pattern of frequency depen¬ 

dence for VLg at small frequencies. This is consistent with our discussion of this model in 
Sect. IV; if the linear frequency dependence sets in at all in this model, that happens only at 
frequencies much smaller than those in Fig. 3. This makes it difficult for us to extrapolate 
our numerical results for this model to the maximal today’s wavelength, /max ~ 10® cm, or, 
equivalently, minimal frequency /min of order 10^ Hz. We have no reason to suppose, how¬ 
ever, that VLgh? can actually increase to small frequencies, beyond its value of order 10“^^ 
at / ~ 10® Hz. 

To conhrm that Fig. 3 is not a numerical artifact, we have made runs in which the 
interaction was switched off, i.e. was set to zero, at r = 100, for the case q = 30. To 
exclude also the effect of self-interaction, the term was replaced at that time by 
The system then becomes a collection of free particles (plus the oscillating zero-momentum 
background) and, if solved exactly in infinite space, should not radiate. In our simulations, 
the intensity of radiation during the interval from r = 100 to r = 125 was three orders of 
magnitude smaller than it was during the same time interval with the interaction present. 


VII. CONCLUSION 

We have shown that gravitational radiation is produced quite efficiently in interactions 
of classical waves created by resonant decay of a coherently oscillating field. For simple 
models of chaotic inflation in which the inflaton interacts with another scalar held, we hnd 
that today’s ratio of energy density in gravitational waves per octave to the critical density 
of the universe can be as large as 10“^^ at the maximal wavelength of order 10® cm. In the 
pure A/^ model, the maximal today’s wavelength of gravitational waves produced by this 
mechanism is of order 10® cm, close to the upper bound of operational LIGO and TIGA 
frequencies. The energy density of waves in this model likely to be well below the sensitivity 
of LIGO or TIGA at such frequencies. 

In other types of inhationary models (or even with other parameters) the ehect can be 
much stronger. We do not exclude that among these there are cases in which it can be 
observable already by LIGO or TIGA. The relevant situations are: 

1. At some values of the coupling constant (or the resonance parameter q), the most 
resonant momenta are close to /c = 0. In the model with massless inhaton this happens, 
for example, for q = 100 (and does not happen for g = 30 or g = 105, which we 
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discussed so far; in these cases, the most resonant momenta are at A; ~ 1). The 
lowest frequency that we had in the box in our numerical simulations for q = 100 
was f ^ 8 X 10® Hz. At that frequency, Qg for q = 100 was almost two orders of 
magnitude larger than for q = 105. In addition, the entire spectrum of gravitational 
waves appears to be shifted to the left with respect to the spectra shown in Fig. |^. 
Cases when resonance is “tuned” to be close to A; = 0 deserve further study. Question 
remains, how large, in such cases, the intensity of gravitational waves can be at the 
horizon scale, Hch- 


2. It is important to consider in detail models of hybrid inflation |^, where the oscillating 
held need not be the inhaton itself, and so the frequency of the oscillations may be 
unrelated to the inhaton mass. 


3. In models where large huctuations produced at preheating cause non-thermal phase 
transitions, as suggested in Ref. [Q], domains or strings can form. A large amount 
of gravitational radiation can be produced in collisions of domain walls, in a way 
somewhat similar to how it happens P,pi]| in models of hrst-order inhation, or in 
decays of a string network, cf. Refs. |M. In particular, in cases when domains are 


formed, intensity of gravitational radiation at the horizon scale, at the moment when 
the domain structure disappears, is expected to be much larger than in cases without 
domains. 
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FIGURES 



T 

FIG. 1. Variances of fields X (solid curve) and (j) (dotted curve) as functions of conformal time 
in the model with massless inflaton for A = 10“^^ and q = 30. 



FIG. 2. Power spectrum of the field (j) for the same model as in Fig. 1, output every period at 
the maxima of (Pq{t)\ k is rescaled comoving momentum (see text). 
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FIG. 3. Today’s spectral density of gravitational waves in the pure A(/)^ model (solid line) and 
in the model where interaction with a massless scalar field X is added; the dashed line 

corresponds to q = 30, and the dotted line corresponds to q = 105, where q = We used 

9*1 m = 100 . 
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